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ABSTRACT 

We discuss the wavelength dependence of angular diameters of M giants from an observational 
perspective. Observers cannot directly measure an optical-depth radius for a star, despite this 
being a common theoretical definition. Instead, they can use an interferometer to measure the 
square of the fringe visibility. We present new plots of the wavelength-dependent centre-to- 
limb variation (CLV) of intensity of the stellar disk as well as visibility for Mira and non-Mira 
M giant models. We use the terms "CLV spectra" and "visibility spectra" for these plots. We 
discuss a model-predicted extreme limb-darkening effect (also called the narrow-bright-core 
effect) in very strong TiO bands which can lead to a misinterpretation of the size of a star in 
these bands. We find no evidence as yet that this effect occurs in real stars. Our CLV spectra 
can explain the similarity in visibilities of R Dor (MSIIIe) that have been observed recently 
despite the use of two different passbands. We compare several observations with models and 
find the models generally under-estimate the observed variation in visibility with wavelength. 
We present CLV and visibility spectra for a model that is applicable to the M supergiant a Ori. 

Key words: techniques: interferometric, stars: fundamental parameters, stars: variables: Miras, 
stars: individual: R Dor, a Ori 



X 



1 INTRODUCTION 

The radius of a star is considered to be one of its fundamen- 
tal properties. Stellar radii may be determined from interfer- 
ometric observations of stellar angular diameters, but stars 
do not necessarily appear to us as sharp-edged disks of uni- 
form brightness and size. Rather, two effects may be present: 
limb-darkening which tapers the surface brightness down to- 
wards the star's edge and line absorption which changes the 
star's apparent size. Both effects depend on the wavelength 
of the observation and are particularly significant in M gi- 
ants. 

To determine the radii of M giants from interferomet- 
ric data, observers have relied upon comparisons with mod- 
els of the centre-to-limb variation (CLV) of the intensity of 
the stellar disk. The simplest CLV, a uniform disk, is often 
used for this. Alternatively we can use CLVs derived from 
appropriate model photospheres of late M giants, since this 
procedure can lead to improvements in these models and a 
greater understanding of these stars. The model-photosphere 
CLVs are defined for a particular filter or spectral passband 
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and must be transformed to visibilities for comparison with 
the observations. Model CLVs for late M giants were first 
studied by Watanabe & Kodaira (1979) and Scholz & Takeda 
(1987; hereafter ST87), and were recently presented for non- 
Miras by Hofmann & Scholz (1998; hereafter HS98) and for 
Miras by Hofmann et al. (1998; hereafter HSW98). 

Model CLVs from these studies have been used for com- 
parison with observations of M giants (for example, Haniff 
et al. (1995), Weigelt et al. (1996) and Perrin et al. (1999)). 
Haniff et al. (1995) have shown that, in some cases at least, 
Gaussian CLVs are better fits to the observations than the 
model CLVs. Other work, e.g. Wilson et al. (1992), has 
shown that some M giants do not appear circular and some 
show asymmetrical brightness distributions. These effects 
have not yet been included in any model. Most observa- 
tions to date have been made at only a few isolated, nar- 
row passbands that were determined by the availability of 
filters, the need to prevent fringe-smearing and the position 
of the deep near-IR absorption features in M giant spectra. 
Recently Burns et al. (1998) and Young et al. (1999) system- 
atically studied the variation in angular diameter with time 
of the Miras R Leo and x Cyg and showed that periodic di- 
ameter variations did occur 



2 A. P. Jacob et al. 



To further improve the M giant models, Scholz (1985) 
and Bessell et al. (1989b) have called for measurements 
of M giant radii over a wide range of wavelengths. Fitting 
wavelength-dependent model CLVs to multi- wavelength in- 
terferometric observations, at a single epoch, should be a 
more robust process for comparing models and observations 
than simply fitting a small number of model CLVs to narrow- 
band observations. This should provide a useful tool for im- 
proving our understanding of M giants, particularly their 
temperature structure. 

In this paper we present the first stage of this work. 
We calculate wavelength-dependent CLVs for several mod- 
els and use them to derive wavelength-dependent visibil- 
ity profiles. We discuss how these visibility profiles can 
be compared with observations. We adopt the perspective 
of an observer who cannot measure an optical-depth ra- 
dius (the usual theoretical radius definition) and has Um- 
ited spatial frequency coverage and therefore limited accu- 
racy in reconstructing the true CLV We discuss an extreme 
hmb-darkening effect that occurs in some models and com- 
pare some past observations with wavelength-dependent vis- 
ibilities derived from appropriate models. As an example, 
we consider the star R Doradus (M8IIIe) which is a semi- 
regular variable with Mira-like properties (Bedding et al. 
1997, 1998). Interferometric observations are usually carried 
out through narrow-band interference filters to avoid fringe 
smearing effects but we have acquired a large set of multi- 
wavelength observations of R Dor using MAPPIT (Masked 
APerture-Plane Interference Telescope, see Bedding et al. 
1994, Jacob et al. 1997) which are currently being processed. 
For these observations MAPPIT included a prism to disperse 
the interference fringes, thus producing multi-wavelength 
interferometric observations. 



2 STELLAR RADH 

Although the stellar radius, R, and its associated effective 
temperature T^s fx (L/i?^)^/"* are convenient and useful 
parameters for describing a star's global properties, limb- 
darkening and Une absorption effects mean we have to 
choose an outer layer whose distance from the centre of the 
star we define to be the radius. Baschek et al. (1991) have 
discussed various radius definitions occurring in the litera- 
ture. The difficulties in formulating meaningful definitions 
are most severe for stars whose outer portions have low den- 
sity gradients and are strongly extended, e.g., Wolf-Rayet 
stars and Mira variables. These stars have extended atmo- 
spheres whose geometrical thickness is not small compared 
to the total dimensions of the stars, in contrast to the compact 
or plane-stratified atmospheres of the vast majority of stars 
in which almost all photons seen by an observer originate at 
essentially the same distance from the star's centre. But even 
main sequence stars such as the Sun present ambiguities, as 
shown by the difference between the photospheric, coronal 
and radio sizes of the solar disk. 

Only for the Sun can we easily measure the brightness 



distribution of the disk. For other stars we can use an in- 
terferometer to measure the (square of the) fringe visibility 
generated by the brightness distribution. Since the visibility 
is usually only measured out to, or just past, the first null and 
with errors of several percent, one cannot uniquely recon- 
struct the CLV. In some cases a null is not present and we are 
limited simply by the available spatial frequency coverage. 
One can, at best, exclude certain extreme CLVs if the data 
are accurate enough (see references given in Scholz 1997; 
hereafter S97, and Hajian et al. 1998). Instead of recon- 
structing the CLV, observations are usually fitted with model 
curves or, if no suitable model is available, with a standard 
limb-darkening profile such as a uniform disk (UD), fully- 
darkened disk (FDD) or a Gaussian. A point somewhere 
along the CLV can then be chosen to be the radius. For the 
Sun the inflection point is used for defining the photospheric 
radius. Thus, the radius we assign to a star is dependent on 
both the model and the chosen radius definition. 



The most common monochromatic radius definition to 

be found in theoretical work is the position of the t\ = 1 
layer, i.e. the radius on a line of sight to the disk's cen- 
tre at which the optical depth reaches unity for a very nar- 
row passband with central wavelength A. However, this is 
not an observable quantity, even for the Sun, and there is 
no simple rule that gives the t\ = 1 point on the CLV. It 
must be be taken from a model. This optical-depth radius 
describes the "typical" layer from which photons seen by 
the observer are emitted. One must be aware that the in- 
tensity contribution function covers a wide range of depths 
around this layer, which may lead, at different wavelengths, 
to quite different CLVs and quite different t\ = 1 positions 
on these CLVs. We refer to S97, HS98 and HSW98 for de- 
tails. For broad-band or whole-spectrum definitions of the 
radius, some method of averaging the optical depth radius 
over wavelength is required. It is common in this case to use 
the Rosseland radius, i.e. the radius at which the Rosseland 
optical depth equals unity (tross = 1). as this is often used 
in modelUng studies. 



The CLV takes into account both Umb-darkening and 
size variation resulting from (molecular) line absorption ef- 
fects. The difficulty in assigning a radius to M giants, both 
Miras and non-Miras, is increased by their strongly extended 
atmospheres and by the almost complete failure of the sim- 
ple UD approximation. Instead, model CLVs similar to FDD 
and Gaussian shapes are more appropriate. In these cases it 
is usual to define a radius to avoid using the fully tabulated 
description of the surface brightness distribution for measur- 
ing the star's "size". On the other hand, since M giants do not 
have a unique radius it may be more useful to compare the 
observed visibility profile with the model visibility profile 
derived from the model CLV. To do this we need to choose 
an appropriate model. 
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3 MODEL CLVS AND VISIBILITIES 

3.1 Models 

We selected from the grid of Bessell et al. (1989a, 1991; 
hereafter collectively BBSW) the IM© static (i.e. non-Mira) 
AGE model, referred to as "X280" in BBSW. This model's 
effective temperature and atmospheric extension (measured 
as the difference between the uniform disk radii in TiO and 
pseudo-continuum bands, as given in HS98) are close to the 
observed values of R Dor. Bedding et al. (1997) found the ef- 
fective temperature of that star to be 2740 K and Jacob et al. 
(1997) estimated the extension to be 20%. Here we use the 
term "pseudo-continuum" to refer to the high-flux spectral 
regions between the TiO bands that lie below about 1 /um, 
since the flux here is still well below the true continuum. The 
BBSW models are static but are suggested by HS98 to be fair 
approximations to semi-regular variables like R Dor, whose 
atmospheric structures are less dominated by outflow and in- 
flow dynamics than those of Mira variables. However, since 
R Dor also shows some Mira-like properties, a IM© Mira 
model, referred to as "E8380" in Bessell et al. (1996; here- 
after BSW96), was also investigated. A further Mira model 
of the P series ("P74200") of HSW98 is used for discussing 
extreme forms of limb-darkening in strong TiO bands. We 
also investigate a model of the H series ("H350") of BBSW, 
which is applicable to the well-studied M supergiant a Ori. 
See Table 1 for basic parameters of the models used. Note 
that although none of the models are tailored specifically for 
the stars listed in the last coluimi they are representative and, 
in some cases, have been used in previous studies for these 
stars. The "parent-star" concept for Mira models is explained 
in BSW96. 

We refer to BBSW, BSW96 and HSW98 for an exten- 
sive discussion of these static and Mira models and the cur- 
rent approximations used in the treatment of molecular band 
absorption, which result in uncertainties in the calculated 
spectrum and CLV in strongly saturated TiO bands. Nev- 
ertheless these models are the most sophisticated presently 
available for these types of stars. 

3.2 Model-Derived CLVs 

For each model we produced CLVs for wavelengths from 
550 to 1100 nm at 1 nm intervals. Each CLV was calculated 
for a rectangular passband of 1 nm width. In order to define 
a reference radius for each of these passbands we carried 
out an intensity-weighted integration of the ta = 1 radii 
over 1 nm in wavelength using the filter radius (TfiUor = 1) 
definition proposed by ST87. For the purposes of this paper, 
we will refer to these narrow passbands, their CLVs and filter 
radii as monochromatic. 

Fig. 1 shows monochromatic CLVs for the IM0 static 
model (X280) at two wavelengths, with normalised intensity 
plotted against radius. Although the intensity scale is nor- 
malised to 1.0 the ratio of the central intensities at pseudo- 
continuum and TiO bands can be very high. In this case 
it is about 2.5, but in some models at some neighbouring 
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Figure 1. Centre-to-limb variations (CLVs) for the IM© static M giant 
model, X280, of BBSW. This shows the effect of both hmb-darkening 
and molecular absorption in a TiO (850 nm) band and a pseudo-continuum 
(820 nm) band as a function of radius (normaUsed to the Rosseland radius). 
The Tfiiter = 1 radius for each passband is marked by a cross. 



TiO/pseudo-continuum bands it can reach 1000. The radius 
scale is normalised to the Rosseland radius for consistency 
with HS98 and HSW98. This normalisation is indepen- 
dent of wavelength. Fig. 1 shows a pseudo-continuum CLV 
at 820 rmi and a neighbouring TiO band CLV at 850 nm. 
The 820 nm CLV is more strongly limb-darkened than the 
850 nm one. In each case the cross marks the rsiter = 1 point 
which, like the t\ = 1 point, is not directly observable. 

Based on this plot we can conclude, from both the over- 
all CLV width and the ta = 1 points, that the star is larger 
in the TiO band than in the pseudo-continuum band. 

In order to show wavelength-dependent CLVs most ef- 
fectively, we have produced CLV "spectra" (i.e. plots of CLV 
versus wavelength). To understand these plots imagine plac- 
ing many monochromatic CLVs (such as those in Fig. 1) be- 
side each other as in Fig. 2. We then choose a fraction of 
the the maximum (i.e. central) intensity, say 50%, (point 1 
on Fig. 2) and find the points on the CLVs corresponding to 
this (point 2). We then project the resulting points onto the 
radius-wavelength plane (point 3) and draw a curve through 
them (point 4). For the case of the IM© static model (X280) 
this produces the line labelled R50 in the top panel of Fig. 3. 
Similarly we can form CLV spectra at 2, 10 and 90 % of the 
maximum intensity (labelled R2, RIO and R90 in Fig. 3). 
This type of plot shows the overall shape and any low-level 
extension of each monochromatic CLV, and also the vari- 
ation in these features with wavelength. The conventional 
intensity spectrum is included in the bottom panel of Fig. 3 
to indicate the position of the various absorption bands. The 
middle panel is discussed in Section 3.3. 
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Table 1. Basic parameters of models: visibility phase, period P, mass M, luminosity L, Rosseland radius of the non-pulsating parent-star of the Mira model 
series Rp, Rosseland radius of non-Mira and Mira model stars R and effective temperature T^g oc (L/R^^^ _^))^''*- The X and H models are static, the 
P and D models are fimdamental mode Miras and the E model is a first-overtone Mira. The last column lists the stars to which the models are compared in this 
paper. 



Model Vis. phase P/days M/Mq L/Lq Rp/Re R^^^^^^^i^ / Rq ^(tb.oss=i) / '^cff /K Reference Star 



X280 








10* 
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a Ori 
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1.00 


328 




6750 


366 




1.09 


2620 


BSW96 


R Dor, R Cas 
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Figure 2. The construction of the CLV "spectrum". See text. 



The shape of a CLV, such as that in Fig. 1, shows how 
the intensity varies as an observer's line-of-sight shifts from 
the disk's centre (where the observer sees deeper layers) to- 
wards the limb (higher layers). This variation includes the 
effects of both limb-darkening and molecular absorption. 
Therefore the R90 and R50 curves in Fig. 3 describe the 
compactness of the CLV. For this slightly-extended model 
atmosphere we see that, above about 750 nm, the brightness 
is more concentrated towards the disk's centre at pseudo- 
continuum wavelengths than in TiO bands. The opposite is 
true below about 750 nm. Only at very low intensity levels, 
say at R2, does one obtain meaningful intensity radii (cf. 
Baschek et al. 1991), i.e. intensity radii that are physically 
related to the typical formation depths of observed spectral 
features. These show that the star is slightly larger in the 
light of TiO photons, that originate in its outer layers, than 
in the light of continuum photons that are formed deeper in 
its atmosphere. The width of the central portions of the CLV 
in either molecular bands or pseudo-continuum wavelengths 



depends on temperatures and temperature gradients and on 
the resulting intensity contribution functions in the relevant 
layers (S97). The assumption that an M giant star always 
appears "larger" when observed in a molecular band is too 
simple. This is discussed further in Section 5. 



3.3 Model-Derived Visibilities 

An observer does not directly see the surface brightness dis- 
tribution of a star, but can measure fringe visibility with an 
interferometer. If we could measure both the amplitude and 
phase of the visibility, and deconvolve the instrumental pro- 
file, we could reconstruct the surface brightness distribution 
from the sampled interferometric data. The phase informa- 
tion tells us about asymmetries in the brightness distribution 
but the only models available are radially symmetric and do 
not include the effects of, for example, hotspots. Therefore 
we will only consider the amplitude of the visibility in this 
paper. 

We converted each monochromatic CLV to a monochro- 
matic visibility curve using the Hankel transform, which 
gives the response of a one-dimensional interferometer to 
a circularly symmetric disk (see Hestroffer 1997). To fix the 
spatial frequency scale we must assume an angular size for 
the model. Following BBSW and BSW96 we have used the 
tross = 1 and the parent-star tross = 1 points (for static 
and Mira models respectively) to define the radius. We as- 
sumed this radius to be 25 mas. This is the approximate size 
of the stars referred to in this paper. Note that the Rosseland 
radius of a Mira model at a given phase is often close, but not 
identical, to the Rosseland radius of the parent-star model. 

Fig. 4 shows the visibility curves for the two monochro- 
matic CLVs of Fig. 1, with visibility plotted against spatial 
frequency. Spatial frequencies are measured in rad~^ and 
they are numerically equivalent, for a 2-element interferom- 
eter, to the baseline length measured in wavelengths. Based 
on this plot we see that the star is larger in the TiO band be- 
cause its visibilities are lower at all spatial frequencies out to 
the first null. This concurs with our conclusion from Fig. 1. 
If we had the complete visibility curves we could also recon- 
struct the true shape of the CLVs. 
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Figure 3. Top: CLV spectrum for the IMq static model, X280. This shows the CLV radius at 2, 10, 50 and 90% of the central intensity, labelled R2, RIO, 
R50 and R90 respectively. Each curve spans 550 monochromatic CLVs. The radius scale is normalised to the tjioss = 1 radius. Middle: Visibility spectra for 
the model at spatial frequencies of 1, 2, 3 and 4x10® rad~^. The Rosseland radius is now assumed to be 25 mas to fix the spatial frequency scale. Bottom: 
Intensity spectrum for the model. The spectrum is plotted twice, with vertical scales differing by a factor of 20, to accommodate the large range in intensities. 
To guide the eye, dashed lines indicate the strongest absorption bands. 
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Spatiol Frequency [l/rad] 

Figure 4. Two visibility curves for the IMq static M giant model, X280, of 
BBSW. These are Hankel transforms of the CLVs in Fig. 1. To fix the spatial 
frequency scale we have assumed that the star has an angular diameter of 
50 mas at the Rosseland radius point in Fig. 1. 

We next used the monochromatic visibility curves to 
produce a visibihty "spectrum" (i.e. visibihty versus wave- 
length) at several different spatial frequencies just as we 
made the CLV spectrum. This is shown in the middle panel 
of Fig. 3. These curves look somewhat like inverted versions 
of the CLV spectra. 

3.4 Broad-band Visibilities 

In order to compare the predictions of the models with ob- 
servations we must include the effect of a real, broad-band, 
filter on the shape of the CLV. This is usually done by in- 
cluding the filter profile in the calculation of each narrow- 
band CLV. However, to make wavelength-dependent CLVs 
for many different filter profiles (or passbands) in this way 
and subsequently make a CLV or visibility spectrum would 
be time consuming. Instead, we produced the monochro- 
matic CLV and visibihty spectra, such as in Fig. 3, and then 
simply smoothed them by convolving each with the required 
observing passband and then normaUsing by the intensity 
spectrum, i.e.. 



where Sconv is the smoothed CLV or visibility spectrum, S 
is the monochromatic CLV or visibility spectrum, / is the 
intensity spectrum, / is the observing passband function and 
(g) denotes convolution. This was done individually for each 
radial position (i.e. for each Rx in Fig. 3) and for each spatial 
frequency. 

As an example we present Fig. 5. Here the IMq static 
(X280) model's CLV and visibihty spectra from Fig. 3 have 



been smoothed, using Eqn. 1, to simulate the effect of a 

4 

super-Gaussian (i.e. ) filter with a 20 nm FWHM. The 
intensity spectrum was smoothed by a simple convolution 
with the same filter profile. This profile approximates quite 
well the interference filters that are often used. The superim- 
posed points from Jacob et al. (1997) are discussed in Sec- 
tion 6.2. 

To check the accuracy of this method we compared 
broad-band CLV and visibility spectra for the X280 model, 
that were produced by smoothing with Eqn. 1, to spectra 
produced directly by including the filter profile in the cal- 
culation of each narrow-band CLV. The difference between 
the CLV spectra produced by these methods was generally 
below 2% and the difference between the visibihty spectra 
was below 0.5%. In Fig. 5 a CLV spectrum at 50% of the 
central intensity, that was produced by the direct method, is 
plotted. It is almost indistinguishable from the R50 curve. A 
visibility spectrum, produced by the direct method, is plotted 
at 2x10^ rad~^ and is indistinguishable from the spectrum 
smoothed using Eqn. 1 . A check using the E8380 model pro- 
duced similar results. 

In going to a smoothed plot. Fig. 5 shows that there 
remain obvious variations in the visibility which should be 
observable. Also, multi-wavelength visibilities as measured 
by MAPPIT provide many more data points for comparing 
models with observations than do present narrow-band tech- 
niques. Clearly this provides a test of how well the models 
predict the variation in visibility with wavelength at a single 
epoch. The narrow-band work done to date, by us and other 
authors, has only been a test of how well the models pre- 
dict the variation of visibility with spatial frequency in that 
band. Although several bands are normally used in narrow- 
band work they are sometimes observed at slightly different 
epochs. In the intervening time the phase, and hence the vis- 
ibility, may change (see e.g. the "D" models in Fig. 11). 



4 THE OBSERVATIONAL PERSPECTIVE 

An interferometer measures the fringe visibility, but only 
to a Umited spatial frequency. This limited range, com- 
bined with measurement errors and the practical difficulties 
involved in detecting high-spatial-frequency low-contrast 
fringes, makes it difficult to directly recreate the true CLV of 
a star Instead observers rely on model fitting. Standard limb- 
darkening profiles (such as UD, FDD or Gaussian) or CLVs 
derived from model photospheres, such as those of BBSW, 
are converted to visibilities and fitted to the observed visibil- 
ity curves. Only at this point can we determine a value for 
the stellar "radius". This could be the UD or FDD radius, 
some reasonable multiple of the HWHM of the Gaussian, 
the T\ = 1, Tfiiter = 1 or tross = 1 optical-dcpth radius, or 
an intensity radius deduced from the shape of the CLV. 

There are several problems that arise in this process. 
Firstly, there are known modelling limitations. Observations 
have shown that M giant and supergiant disks exhibit bright 
surface features and/or geometric elongations (see e.g. Wil- 
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Figure 5. As for Fig. 3 but smoothed to simulate tlie effect of observing with a 20 nm I^WHM super-Gaussian passband. The dotted curve overlying the R50 
CLV spectrum is discussed in Section 3.4. Superimposed on the visibiUty spectra are data for R Dor (asterisks and diamonds) interpolated from Fig. 3 of Jacob 
et al. (1997). These data are discussed in Section 6.2. 



son et al. 1992), but no model yet incorporates these phe- 
nomena. The BBSW models also have limited reliability in 
deep absorption bands. 

Secondly, most of the BBSW and BSW96 model CLVs 



show a null in their visibilities at spatial frequencies acces- 
sible to 4m-class telescopes but most observations do not 
show it (see e.g. Haniff et al. 1995, Weigelt et al. 1996, Bed- 
ding et al. 1997). Instead, Gaussian curves provide better fits 
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to the observations. Haniff et al. (1995) suggested that this 
could be due to localized hotspots or circumstellar scattering 
that increases the visibility. If this is the case then comparing 
multi-wavelength visibilities with model predictions, at low 
spatial frequencies which are relatively unaffected by these 
phenomena, could provide a way to improve the models un- 
til these effects are included. 

Thirdly, an observer's interpretation of the data is hin- 
dered by the similarity of the visibility curves within the first 
null even when the CLVs are significantly different. For ex- 
ample, in Fig. 3 the visibilities between 1000 and 1100 nm, 
which includes a VO band, are approximately constant with 
wavelength, even though the CLV shape varies significantly 
across this band. In this case it is the limb-darkening that dif- 
fers. Fig. 6 shows the CLV, visibility and intensity spectra for 
a IM0 Mira model. Here the visibility spectrum has a dip at 
the VO band. Since the models are expected to give reason- 
able predictions at these wavelengths (Brett 1990), this part 
of the visibility spectrum could be useful in distinguishing 
between different models, despite the difficulty in observing 
these wavelengths. Additionally, the similarity in the visi- 
bility curves means that if the diameter is treated as a free 
parameter it can be impossible to tell whether two visibil- 
ity curves are different due to different diameters (given a 
single limb-darkening profile) or to different limb-darkening 
profiles themselves, unless one looks beyond the first null. 

Finally there is the possible misinterpretation of model 
visibilities by an observer. A low visibility within the first 
nuU, at a particular spatial frequency, means the interferom- 
eter is resolving the source more than if the visibility was 
high. We interpret this to mean the source is larger in angu- 
lar size. If we assume the TiO overlies the pseudo-continuum 
forming layers, as it must do to form absorption bands, then 
we would expect to see a larger star in the TiO bands than 
in the pseudo-continuum bands. Observationally we would 
expect lower visibilities (within the first null) in TiO bands. 
However, in some models in some TiO bands we noticed 
that the visibility was higher than in neighbouring pseudo- 
continuum bands. See, for example, the visibilities below 
about 750 nm in Fig. 3 and Fig. 6. A first glance at the model 
visibility curves might lead an observer to conclude that the 
star is smaller in these TiO bands than in the neighbouring 
pseudo-continuumbands. This does not seem to make sense. 
The explanation Ues in the form and extent of the limb- 
darkening present (see section 5). The observer's misinter- 
pretation of these model visibilities results from not consid- 
ering carefully the visibilities near and beyond the first nuU 
or at high spatial frequencies. 

Clearly the visibiUty spectrum, and monochromatic vis- 
ibiUty curves, require careful interpretation. 



5 AN EXTREME LIMB -DARKENING EFFECT IN 
DEEP TITANIUM-OXIDE BANDS 

A star with a compact atmosphere, such as the Sun, has a 
relatively sharp edge to its brightness distribution. The dis- 



tance of this edge from the disk's centre is effectively inde- 
pendent of wavelength and thus defines the stellar radius. If 
observations of a compact-atmosphere star appear to imply 
wavelength-dependent radii this indicates incorrect assump- 
tions have been made about Umb-darkening. In fact, this pro- 
vides a test for whether or not we have chosen the correct 
limb-darkening form at each wavelength. For example, if a 
single limb-darkening model (e.g. UD) is fitted to observed 
visibilities at two wavelengths that actually have different 
CLVs this will, incorrectly, give different angular diameters. 

In the case of an extended-atmosphere star this effect 
may occur in conjunction with a genuine physical differ- 
ence in the extension of the stellar disk seen at different 
wavelengths. The CLV in a strong-absorption spectral fea- 
ture formed in high atmospheric layers (i.e. large t\ = \ 
radius) may be flatter than the CLV in a pseudo-continuum 
feature formed in deep layers (small t\ = 1 radius). For a 
small atmospheric extension, such as occurs in most static M 
giants, the CLV behaviour (e.g. Watanabe & Kodaira 1979, 
ST87, HS98) may obscure the larger size of the star at the 
strong-absorption wavelength. Nevertheless, the geometric 
extension effect normally prevails in very cool non-Miras 
and in Miras having substantially extended atmospheres. In 
this case the central maximum of the visibility is narrower at 
the wavelength of the strong-absorption feature which indi- 
cates a larger disk. Data reduction by means of an identical 
limb-darkening model would lead to an underestimate of the 
ratio of the disk sizes in strong and weak absorption, but this 
would not suggest a smaller disk at the stronger spectral ab- 
sorption feature. 

We found, however, during the present investigation that 
some of the BSW96 and HSW98 Mira models predict ex- 
treme CLV shapes in very strong TiO bands. We use the 
term "extreme" here for the "narrow-bright-core" effect de- 
scribed in HSW98. To save space not all of these models are 
presented here, although the IMq Mira model (E8380) in 
Fig. 6 is a typical example. The width of the central maxi- 
mum of the visibility in very strong TiO bands (e.g. 710nm) 
becomes significantly broader than that seen at neighbour- 
ing pseudo-continuum wavelengths, a situation that at first 
glance seems incorrect. The disk, if seen directly by an ob- 
server, would be extremely limb-darkened and have very ex- 
tended faint outer portions. The CLV has a Gaussian shape 
and the t\ = 1 optical-depth radius is positioned in the low- 
intensity wing (see e.g. Fig. 8 of HSW98). Only very high 
accuracy measurement of the visibility curve at high spatial 
frequencies would reveal the large size of the stellar disk in 
the TiO band. 

This extreme limb-darkening effect may also show up 
in some very cool non-Mira models (see Fig. 3). However, 
it is most pronounced in certain Mira models such as the 
near-maximum IM© model (E8380). Fig. 6 shows the ex- 
treme central concentration of the strong-TiO band bright- 
ness distributions that occur below about 750 nm in this 
model. Other extreme examples are found at near-minimum 
phases of the P series of models of HSW98 (see the P71800 
model in Fig. 2, upper left panel, of S97). However, the near- 
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Figure 6. As for Fig. 3 but for the IMq Mira model, E8380. Here the radius is normalised to the parent star's Rosseland radius, Rp. 



maximum model (P74200) of the same series (Fig. 2, upper 
right panel, of S97) does not show this effect. Fig. 7 shows, 
for this near-maximum model, the "expected" situation, i.e. 
the central visibility maxima in strong TiO bands are nar- 
rower than those in neighbouring pseudo-continuum bands. 



Note, all Mira model CLVs in this paper are displayed with 
their radii normalised to their parent star's Rosseland radius. 

Unfortunately, this effect depends not only on the spe- 
cific model and model phase but also on subtle details of the 
structure of the upper layers (HSW98). Present Mira mod- 
els are still approximate. The temperature stratification in 
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Figure 7. As for Fig. 3 but for the IMq Mira model, P74200. Here the radius is normalised to the parent star's Rosseland radius, Rp. The extreme limb- 
darkening effect is not present but molecular absorption has caused large differences in the star's size between absorption and pseudo-continuum bands. 



the upper atmosphere is significantly affected by the simple 
treatment of TiO molecular band absorption which is par- 
ticularly sensitive for strongly saturated bands (see BBSW, 
BSW96). This simple treatment also does not account for 
the effects of Doppler- shifting in layers with outflowing and 



infaUing matter. Furthermore, the adopted absorption coeffi- 
cients directly affect the CLV shape. 

A search of the literature revealed only two sets of pub- 
lished visibility curves relevant to this discussion. Fig. 2b of 
Haniff et al. (1995) shows visibilities for R Cas at 700nm 
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and 710nm obtained in September 1993. Fig. 4 of Weigelt 
et al. (1996) also shows visibilities for R Cas at 700 nm and 
714nm obtained in September 1994. There is no evidence 
for an extreme limb-darkening effect in this star. We have 
MAPPIT data for R Dor covering these two bands which we 
are currently processing. Therefore, it is not clear at present 
whether this extreme limb-darkening effect does occur in 
real Mira stars or is only a modelUng artifact produced by 
inadequate modelling techniques. 

A demonstration of the influence of the temperature gra- 
dient on the brightness profile is given in Fig. 8. We calcu- 
lated CLVs for the Mira model (E8380) in two rectangular 
passbands of 1 nm width. One was centered at 701.5 nm in a 
pseudo-continuum band and the other at 710nm in a very 
strong TiO band. Fig. 8a shows that the very strong TiO 
CLV (dashed curve 1) is much narrower than the pseudo- 
continuum CLV (solid curve). In Fig. 8b the visibility at the 
very strong TiO band (dashed curve 1) is broader than that 
at the pseudo-continuum band (sohd curve). This is what we 
would expect after seeing Fig. 8a, yet an observer, who only 
sees the visibility curves and usually only within the first 
nuU, might naively conclude that the star was larger in the 
pseudo-continuum band. 

The respective t\ = 1 radii are 2.29 and 1.56 xRp 
(marked by crosses) with the former lying in the very dim 
outer portions of the CLV curve. So, by the ta = 1 def- 
inition, the TiO band radius is larger than the pseudo- 
continuum radius. Note that the ri.o4 = 1 continuum radius 
is at 1.08 xRp (HSW98). 

We now perform a simple numerical experiment to in- 
vestigate the effect on the brightness profile of changing the 
temperature gradient in the outer layers. We assume that the 
outermost atmosphere is optically thin with respect to fea- 
tures affecting radiative equilibrium above 2.0, 1.9 and 1.85 
xRp (where T = 1090, 1180 and 1250 K) sequentially. This 
has the effect that the temperature gradient becomes flatter 
and follows a T{r) oc r~^/^ law. We may consider this as, 
for example, a very simplistic correction for possibly exces- 
sive near-surface cooling which could be caused by over- 
blanketing due to inadequate molecular band treatment. The 
result is shown in Fig. 8. The 701.5 nm pseudo-continuum 
CLV (FDD-like shape) is not affected by these manipula- 
tions. On the other hand, the width of the very-strong-TiO 
CLV (Gaussian shape) becomes progressively broader with 
a corresponding change in the visibility curve (dashed curves 
2 to 4). Flattening the gradient above any layer that lies in- 
side about 1.90 xi?p (dashed curve 2) results in the "ex- 
pected" situation where the strong-TiO visibility is narrower 
than the pseudo-continuum visibility. 

It is obvious from the above that accurate interferomet- 
ric observations of strong TiO bands would provide a sensi- 
tive method of probing the temperature structure of the out- 
ermost atmospheric layers. 



6 COMPARISON WITH PUBLISHED 
VISIBILITIES 

6.1 The 850 nm TiO Band in R Dor 

Bedding et al. (1997) measured visibility curves for R Dor 
using an interference filter of 40 nm width centered at 
855 nm. The pulsation phase of the star at that time is dif- 
ficult to determine due to the recent irregularity of R Dor 
(see Bedding et al. 1998). Jacob et al. (1997) found very 
similar visibilities to the Bedding et al. (1997) results but 
at a pulsation phase of about 0.7 and with a bandwidth of 
20 nm centred at 850 nm. These visibilities agreed very well 
with each other despite the different observing passbands. 
Both the static (X) and Mira (E) models used in this paper 
to represent R Dor can be used to explain this agreement. 
Fig. 3 and Fig. 6 show that the visibility spectra at all spa- 
tial frequencies are approximately constant from 820 rmi to 
870 nm. This is also the case for other phases of the E model. 
Thus, as long as the wings of the filter are within this region, 
the measured visibiUty will be independent of the filter width 
and central wavelength. 

In fact the visibility spectrum could be a useful tool for 
choosing the filter width, its central wavelength and the al- 
lowable leakage from, for instance, an absorption band. Ob- 
servers normally use the intensity spectrum to choose filter 
properties. 

6.2 Otlier Multi-wavelengtli Visibilities 

We are aware of only four sets of published multi- 
wavelength visibilities that we can compare with the BBSW 
models using our visibiUty spectra. They are for R Dor (Ja- 
cob et al. 1997), R Cas (Haniff et al. 1995 and Weigelt et al. 
1996) and R Leo (TuthiU et al. 1994). 

Our criterion for choosing these sets of observations 
was that they include adjacent TiO/pseudo-continuum bands 
at the same epoch (i.e. same pulsation phase). Figs. 5 and 
9 to 12 show these observations over-plotted onto rele- 
vant models. In each plot, the visibility spectra have been 
smoothed, using Eqn. 1, to closely match the FWHM of the 
filters used for the observations. The data points are Unearly 
interpolated from the relevant figures in each paper where 
necessary. We have not included errors in these plots but the 
errors in the measured data ranged from 6 to 20 %. Straight 
lines between data points simply connect points at the same 
spatial frequency for clarity. The model visibilities are plot- 
ted for either a Rosseland (for static models) or parent-star 
Rosseland (for Mira models) radius of 25 mas for simplicity, 
i.e. no attempt has been made to match the size of the model 
star to the size of the observed star. Therefore, this compari- 
son shows to what degree the models predict the variation in 
visibility with wavelength, but not how well they predict the 
absolute visibility. 

The Jacob et al. (1997) data for R Dor are plotted on the 
visibility spectra for the IM© static model (X280) in Fig.5. 
The asterisks are observations at 820/18 nm and the dia- 
monds are observations at 850/20 nm. The upper and lower 
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Figure 8. CLVs (a) and visibilities (b) for the IMq Mira model, E8380, at 701.5 (pseudo-continuum) and 710nm (TiO) showing how the TiO band CLV and 
visibility change as the temperature gradient in the outer layers is artificially flattened (see text). Crosses mark the ta = 1 radii points on the CLVs and the 
n.04 = 1 continuum radius is marked by a short dash at 1.08 x Rp. 



pairs of points are at spatial frequencies of 3 and 4 x 10^ 
rad~^ respectively. The model has a much smaller vari- 
ation in visibiUty with wavelength than the observations. 
Since R Dor shows some Mira-like properties we com- 
pare the same data to the IM© Mira model (E8380) of 
BBSW in Fig. 9. Only data points at a spatial frequency of 
3 X 10^ rad~^ have been plotted here and the model has 
been smoothed to simulate the effect of a 20 nm FWHM 
super-Gaussian filter. Again the model has a much smaller 
variation in visibility with wavelength than the observations. 

Haniff et al. (1995) observed the Mira variable R Cas 
at 700/10 and 710/lOnm in September 1993 at a phase of 
about 0.03. In Fig. 10 we compare these observations to a 
fundamental-mode Mira model (D28760) at phase 0.0 (top) 
and to a first-overtone mode Mira model (E8380) at phase 
0.0 (bottom). These models have a phase close to the obser- 
vations and were used in Haniff et al. (1995) to represent 
R Cas. The upper and lower pairs of points in each panel are 
at spatial frequencies of 2 and 3 x 10^ rad~^ respectively. 
Again the observed visibilities vary significantly more than 
the model visibilities. In fact, the models predict the reverse 
of the observed variation. 

Weigelt et al. (1996) observed R Cas at 700/6 and 
714/6 nm in September 1994 at a phase of about 0.91. In 
Fig. 11 we compare these observations to fundamental-mode 
Mira models at phase 1.8 (D28320: top) and phase 0.0 
(D28760; middle) and to a first-overtone mode Mira model 
(E8380) at phase 0.0 (bottom). Note that both Mira mod- 
els and Miras themselves show cycle to cycle variations in 
their behaviour hence pulsation phases can be greater than 
1 .0 (BSW96, HSW98). For the two fundamental-mode mod- 
els the top, middle and bottom pairs of points in each panel 



are at spatial frequencies of 2, 3 and 4 x 10® rad^^ respec- 
tively. For the first-overtone model only the observations at 
spatial frequencies of 2 and 3 x 10^ rad~^ are shown. Only 
the fundamental-mode model approaching maximum bright- 
ness in its second cycle (i.e. at phase 1.8) has a similar vari- 
ation in visibiUty with wavelength to the observations. The 
fundamental model at maximum brightness (i.e. phase 0.0) 
and the first overtone model are not well matched to the ob- 
servations and again show the reversed behaviour. 

Tuthill et al. (1994) observed the Mira variable R Leo 
at 833/41 and 902/50nm in January 1992 at a phase of 
about 0.26. In Fig. 12 we compare these observations to a 
fundamental-mode Mira model (D28847) at phase 0.2 (top) 
and to a first-overtone mode Mira model (E8560) at phase 
0.21 (bottom). We have smoothed the visibility spectra here, 
using Eqn. 1, assuming an intermediate FWHM of 45 nm. 
The upper and lower pairs of points in each panel are at spa- 
tial frequencies of 2 and 3 x 10^ rad~^ respectively. Both 
model and observations show similar variations in visibility 
with wavelength at 2 x 10^ rad~^, although it is difficult to 
say more because of the large FWHM, the relative flatness 
of the visibility spectra here and the position of the 833 nm 
filter being just outside the centre of the TiO band. 

In summary, only one model shows a similar variation in 
visibility with wavelength to the observations. In most cases 
the models underestimate this variation and in several cases 
it is in the opposite sense to the observations, i.e. the TiO 
band (approx. 710nm) visibility is higher than the pseudo- 
continuum band (700 nm) visibility in the models. However, 
this TiO band is subject to the extreme limb-darkening effect 
and the models may be unreliable here. This small sample 
leaves it unclear whether the models accurately predict the 
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Figure 10. Observations of R Cas by Haniff et al. (1995), compared with 
two solar-mass Mira models of different pulsation modes: fundamental 
(top) and first-overtone (bottom). 



variation in visibility with wavelength for M giant stars. Our 
multiwavelength data for R Dor and other stars will provide 
a better test of the accuracy of these models. 



7 AMODELFORBETELGEUSE 

Fig. 13 shows the CLV, visibility and intensity spectra for 

the 15Mq model (H350) which is used to represent the su- 
pergiant a Ori (Betelgeuse, M2 lab). Unfortunately we have 
found no suitable observations (using the criterion in Sec- 
tion 6.2) to compare with this model. 

The model displays much smaller variations in both 
CLV and visibility with wavelength than later M giant mod- 
els. In particular, the R2 radius is independent of wave- 
length, which reflects the fact that the H350 model atmo- 
sphere is nearly compact. Above about 900 nm, its CLV and 
visibility spectra become constant. Below about 800 nm, we 
once again see higher visibilities, within the first null, at 
TiO bands than in pseudo-continuum bands. A naive in- 
terpretation of this, without consideration of the detailed 
CLV shape, would be misleading. Therefore an observer 
of the H350 model star might conclude that the stellar ra- 
dius is smaller at TiO than at pseudo-continuum wave- 
lengths whereas, in reality, the atmosphere is compact and 
diameters are not wavelength-dependent. It is just that the 
limb-darkening effect is very pronounced in the TiO bands. 
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Figure 11. Observations of R Cas by Weigelt et al. (1996), compared with 
three solar-mass Mrra models of different pulsation modes: fundamental 
(top & middle) at different phases and first-overtone (bottom). 



In fact Cheng et al. (1986) explained the scatter in previ- 
ous wavelength-dependent angular diameter measurements 
of a On by adopting a varying Umb-darkening behaviour. 
Although this effect appears similar to the extreme limb- 
darkening effect discussed in Section 5, we prefer to reserve 
that term for cool giants, in particular Miras, that have non- 
compact atmospheres. As discussed in that section, limb- 
darkening in the strongest TiO bands below 800 nm for some 
of those stars may be so extreme that it increases the visibil- 
ities strongly within the first null and disguises the effects of 
geometric extension of the star's atmosphere. This peculiar 
situation is called the extreme limb-darkening effect in this 
paper. 
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Figure 12. Observations of R Leo by Tuthill et al. (1994), compared with 
two solar-mass Mira models of different pulsation modes: fundamental 
(top) and first-overtone (bottom). 



8 CONCLUSION 

Interferometric observations of M giants give us visibility 
profiles related to their centre-to-limb variations (CLVs). 
Observers never see the CLV directly and cannot measure an 
optical-depth radius (which is commonly used in theoretical 
work). They can only measure the visibility with a limited 
accuracy and to a limited spatial frequency which prevents 
reconstruction of the true CLV. Therefore observers convert 



standard limb-darkening profiles or CLVs from model pho- 
tospheres of M giants to model visibihty profiles and fit these 
to the observed visibility profiles. Only then can they de- 
fine a radius, effective temperature and other parameters for 
the star. The MAPPIT interferometer is capable of record- 
ing multi-wavelength visibilities. Comparing these with the 
models should help improve the models and our understand- 
ing of M giants, particularly their temperature structure. 

We have developed a way of displaying wavelength- 
dependent CLVs and visibilities. Our CLV and visibility 
spectra, as we have called them, show CLV versus wave- 
length for several radial positions and visibility versus wave- 
length for several spatial frequencies. The CLV spectra show 
both the overall shape and any low-level extension of the 
CLV at each wavelength while the visibility spectra allow 
direct comparisons with multi-wavelength observations. The 
CLV data are derived from the Mira and non-Mira M gi- 
ant photospheric models of BBSW, BSW96 and HS98. The 
data are displayed at wavelengths from 550 to llOOnm but 
this could easily be extended to cover the range present in 
the models which is 320 nm to 4100nm. We are also able 
to plot other radial positions or spatial frequencies and can 
easily and quickly simulate the effect of any real filter on 
the monochromatic CLV and visibihty spectra. The visibil- 
ity spectra provide a tool for comparing multi-wavelength 
interferometric observations of M giants with models of their 
centre-to-limb variation at a single epoch. We developed 
these plots to help interpret multi-wavelength interferomet- 
ric observations that we have obtained with MAPPIT. The 
plots might also be useful when choosing filter properties 
(central wavelength, width and leakage) for narrow-band ob- 
servations. 

Care must be taken when making comparisons be- 
tween interferometric observations and models. Several ef- 
fects, including the model limitations and an extreme limb- 
darkening effect could lead to misinterpretations of the data. 
The extreme limb-darkening effect occurs in very strong 
TiO bands but only in some Mira models and some very 
cool non-Mira models. To an observer the model visibiUties 




Figure 13. As for Fig. 3 but for the 15M0 static M giant model, H350. The flatness of the R2 and RIO CLV spectra is a consequence of the compactness of 
this model atmosphere. 



might suggest, at first glance, that the star appears smaller 
in the very strong absorption bands. Careful examination 
of high spatial frequency model-derived visibilities, or the 
model CLV, would show the star actually has a faint ex- 
tended Umb. The concept that an M giant wiU always ap- 



pear larger, to an observer, in an absorption band may be too 
simple, at least when examining models. However, the few 
pubhshed observations that are relevant to this situation do 
not show this extreme effect, so further observations would 
be useful in establishing whether or not it is only a mod- 
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elling artifact. Since a star's CLV depends on temperatures 
and temperature gradients, accurate interferometric observa- 
tions of strong TiO bands would provide a sensitive method 
of probing the temperature structure of the outermost atmo- 
spheric layers and improving the models. 

Our visibility spectra can now explain the similarity in 
visibilities of R Dor observed by Bedding et al. (1997) and 
Jacob et al. (1997) despite the different filter properties used 
in those studies. We chose appropriate models from BBSW 
and BSW96 and found that the visibility spectrum is flat near 
850 run and so a range of filter profiles will result in the same 
visibility profile. 

At present there are only a limited number of obser- 
vations that can be compared to models using our visibil- 
ity spectra. In most cases the models show a much smaller 
variation in visibility with wavelength than the observations, 
and in some cases the models predict the reverse of the ob- 
served variation. Only one model, the fundamental-mode 
Mira (D28320), shows a similar variation in visibility with 
wavelength to the observations of R Cas by Weigelt et al. 
(1996). Our multi-wavelength observations of R Dor, which 
are currently being processed, will provide a clearer test of 
how well the models predict the variation in visibility with 
wavelength for M giant stars. 
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